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METHODS AND ARRANGEMENT FOR CHANNEL ESTIMATION IN A 
TELECOMMUNICATIONS SYSTEM 

FIELD OF INVENTION 

The present invention relates to communication, and primarily 
to wireless communication. 

5 DESCRIPTION OF RELATED ART 

In a radio communication system, a bitstreaiti or a symbol stream 
is sent on a radio channel from a transmitter to a receiver. 
Physical limitations of wireless channels, such as bandwith 
limitations, propagation loss, interference, multipath fading, 
10 inter symbol interference and time variance, present a 
fundamental challenge to reliable communication. Further 
challenges come from power, limitation, size, and speed of 
devices that are employed within portable wireless apparatus. 

It is well known for a man skilled in the art that in most 
15 scattering environments, antenna diversity, i.e. the usage of 
multiple transmitting and/or receiving antennas, is the most 
practical technique for reducing the effect of multipath fading. 
The diversity technique involves the provision of a number of 
replicas of the transmitted signal to the receiver, some 
20 replicas being less attenuated by fading 

In a conventional serial data system, the symbols are 
transmitted sequentially on a channel, the frequency spectrum of 
each data symbol being able to occupy the entire available 
bandwith. 

25 A problem with serial transmitting methods in radio 
communication systems is that multipath reception having 
different delayed signals, due to different transmission paths, 
may result in Inter Symbol Interference, ISI. This occurs when 
signals being delayed, due to e.g. reflections, and having a 

30 non-negligible strength, are received having a time difference 
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as long as a data symbol or longer. The most common solution to 
this problem is the usage of an equalizer at the receiver. Based 
on estimates of channel attenuation, delay and phase of a 
channel, the equalizer tries to compensate the received signal 
5 for impact from the channel. In the time domain, the conplex 
attenuation and delay description of a channel is called channel 
impulse response. The channel impulse response also enables the 
measurement of how the received energy is scattered over time by 
calculating a delay power spectrum. The range of time of which 

10 the delay power spectrum is essentially zero is called delay 
spread. A problem in the case of single carrier in a serial 
system, i.e. when the symbols are transmitted one by one after 
each other, is that for a given delay spread, the Inter Symbol 
Interference is provided over several symbols, requiring a more 

15 complex equalizer than in systems wherein the symbols are sent 
concurrently. High complexity in the algorithms requires more 
processing at the receiver, and thus high battery consumption, 
high clock rates, larger chip area, etc... It is therefore of 
great importance to keep the complexity low. 

20 A parallel data transmission system offers possibilities for 
alleviating many of the problems encountered with serial 
systems. A parallel system is one in which several sequential 
streams of data are transmitted simultaneously, so that at any 
instant many data elements are being transmitted. In such a 

25 system, the spectrum of an individual data element normally 
occupies only a small part of the available bandwith. This is 
often denoted multi-carrier transmission. 

Furthermore, the impact of a sudden short lasting impulse 
interference, e.g. sparks, is alleviated by the fact that the 
30 disturbance is spread over multiple channels, as the disturbance 
on each symbol is reduced* Hence, due to the dividing of an 
entire channel bandwith into many narrow sub-bands, the 
frequency response over each individual sub-band is relatively 
flat. Since each sub-channel covers only a small fraction of the 
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original bandwith and hence the channel is flat, an equalization 
method is potentially simpler than in a serial system having 
ISI. A simple equalization algorithm and the implementation of 
differential encoding may make it possible to totally avoid 
5 equalization. 

The invention is directed towards a system using Orthogonal 
frequency division multiplexing OFDM that can be simply defined 
as a form of multi channel modulation where its channel spacing 
is carefully selected so that each sub-channel is orthogonal to 
10 the other sub-channels. Thus parallel channels may be called 
mul t i - c ar r i er . 

In orthogonal frequency division multiplexing (OFDM) , a channel 
is divided into many narrow parallel sub-channels, thereby 
increasing the time of an individual symbol and reducing or 

15 eliminating the inter-symbol interference (ISI) caused by the 
multipath environments. On the other hand, since the dispersive 
property of wireless channels causes frequency selective fading, 
there is a higher error probability for those sub-channels when 
deep fades occurs. Hence, techniques such as error correction 

20 codes and diversity have to be used to compensate for the 
frequency selectivity • 

In order to provide good channel quality estimation, a pilot 
signal, i.e. a training sequence, is used. Here follows a more 
detailed disclosure on how a possible OFDM signal is transmitted 
25 and corrected for channel impairments. 

In this disclosure, capital letters are used to denote frequency 
domain signals whereas lower case letters denote time domain 
signals. Furthermore, the indexes k and n are used for frequency 
and time domain indices, respectively. 

30 First, for the purpose of coherent detection, a pilot signal is 
transmitted. A pilot is a sequence of symbols that are defined 
and well known in both transmitter and receiver. The pilot is' 



WO 02/45329 



4 



PCT/SE01/02292 



denoted P(k) and is a discrete sequence of complex symbols, 
which are indexed in frequency with the parameter k and ranges 
from 0 to N-l. Prior to the transmission, the pilot P(k) is 
transposed into the time domain by an Inverse Discrete Fourier 
5 Transform IDFT resulting in a complex valued discrete pilot 
sequence in the time domain denoted p(n), where n is a time 
index ranging from 0 to N-l. It should be noted that 
practically, an (Inverse) Fast Fourier Transform is used instead 
of an (Inverse) Discrete Fourier Transform, because it has lower 

10 algorithmic complexity and therefore requires less complex 
hardware. The pilot sequence is sent over the wireless channel 
modeled by the discrete and complex valued channel impulse 
response h(n) # which performs a linear filtering function of the 
sequence p(n) . The resulting received sequence is r p (n), which 

15 is discrete and complex valued: 

r p (n)=p(n) *h(n) + v(n) , where 

v(n) is the noise added in the receiver or other ambient noise 
at the antenna input. The filtering process is identified by the 
sign * between the sequence p(n) and the discrete and complex 
20 channel impulse response h(n) . It is a so-called convolution. 

In an OFDM system, the signals are mainly processed in the 
frequency domain. Therefore, the received N sample sequence is 
Discrete Fourier Transformed (the opposite of Inverse Discrete 
Fourier Transform) resulting in 

25 R p (k)=P(k) -H(k)+V(k), where 

the sequence P(k) is a known sequence, whereas V(k) is a random 
sequence and therefore best modeled as a stochastic process. 
H(k) is a channel transfer function which is unknown prior to 
the transmission. V(k) is noise and therefore random. 

30 It should be noted that the convolution * is transformed into a 
multiplication under the Discrete Fourier Transformation. Due to 
implementation reasons (i.e. complexity) a Discrete Fourier 
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Transform is seldom performed but instead an equivalent lower 
complexity, so called Fast Fourier Transform, is carried out. 
The equivalent is done for Inverse Discrete Fourier Transform, 
i.e. the Inverse Fast Fourier Transform is used. 

5 An estimate of the function H(k), which is required for an 
equalization step, is searched. Based on the received signal 
Rp(k), and prior knowledge of the sequence P(k) the function 
H(k) can be determined according to: 

R D (k) 
P(k) 

10 It should however be noted that this is only an estimate, since 
noise is present. If the noise level is low in comparison to the 
received signal of the sequence Fp(k), the estimation of the 
channel is precise. On the other hand, if the noise level is 
high, the estimation is poor. 

15 When the channel has been estimated, data can be transmitted. 
The data sequence (discrete and complex valued) is denoted S(k), 
where k is the frequency index ranging from 0 to N-l. Prior to 
the data transmission, the sequence S(k) is Inverse Discrete 
Fourier Transformed, or preferably Inverse Fast Fourier 

20 Transformed to a time representation s (n) , which is discrete and 
conplex valued. It is assumed that the time elapsed since the 
channel estimation was performed is so short that the channel 
characteristics have not undergone any significant change. 

Accordingly, when the sequence s(n) is transmitted, it will 
25 experience the same channel characteristics as the pilot signal 
did. Of course, the noise is different at reception, which is 
here indicated with w, instead of v. 



r B (n) =s (n) *h (n) +w (n) 
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When the received signal has been Discrete Fourier Transformed, 
or preferably Fast Fourier Transformed, the following signal is 
obtained 

R s (k)= S(k) • H(k) + W(k) 
5 As the channel is known, the sequence S (k) can be estimated as 

;„.s (fc)-^(fc) 

H(k) 

If desired, more data symbols can now be sent- At each time 
instant, the channel estimation is used to compensate for 
channel impairments, i.e. equalization. However, the channel 
10 characteristics are likely to change over time and therefore may 
new pilot signals be sent regularly having a fixed time 
interval . 

However, a man skilled in the art understands that the 
disclosure above is related to an extremely trivial case. In 

15 practice, this is probably not the way equalization would be 
performed. Instead, equalization is often an integral part of 
forwarding error correction decoding. Furthermore, a particular 
a priori knowledge of the channel, such as delay spread or 
statistical characteristics, can be used to improve the quality 

20 of the channel estimate which in turns improves the quality of 
the estimated data S(k) . 

Transmitter diversity, i.e. the usage of at least two 
transmitter means, such as transmitting antennas, is an 
effective technique for combating fading in mobile wireless 
25 communications, especially when receiver diversity, i.e. the 
usage of at least two receiver means, such as antennas, is 
expensive or impractical . 

The performance gain of linear transform based diversity with 
ideal maximum likelihood sequence estimation (MLSE) and an 
30 arbitrary number of transmitter antennas is investigated and 
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compared to receiver diversity in the article by J.H. Winters 
"The diversity gain of transmit diversity in wireless systems 
with Rayleigh fading", Proc. 1994 IEEE Int Communications Conf, 
Chicago, IL, June 1994, pages 1121-1125. 

5 The published International patent application WO 99/14871 
discloses a transmit diversity method for wireless 
communication. In one illustrative embodiment, two transmit 
antennas and a single receive antenna are used, providing, the 
same diversity gain as maximum- ratio receiver combining scheme 
10 with one transmit antenna and two receive antennas. 

More recently, space time coding has been developed for high 
data rate wireless communication. Space time coding schemes are 
based on coding signals in time as well as space, i.e. over 
multiple transmit antennas. Typically, a number M of transmit 

15 antennas and a number N of receive antennas are employed. In the 
article by G. Foschini, "Layered space-time architecture for 
wireless communication in a fading environment when using multi- 
element antennas", Bell Labs Technical Journal Autumn 1996, the 
author describes a multiple transmit and receive antenna system 

20 for which the channel capacity scales linearly with the minimum 
number of employed antenna at the transmit and receive side. It 
should be noted that transmitter diversity is sometimes 
classified under the field of space time coding although merely 
a single antenna is employed at the receiving side. In the 

25 latter context, it is clear that space time coding may 
alternatively to high channel capacity instead be exploited for 
increased robustness performance under e.g. fading conditions as 
transmit diversity methods mainly address. 

In the article by D. Agarwal et. al. * Space- time coded OFDM for 
30 high data rate wireless communication over wideband channels", 
Proc. 48 th IEEE Vehicular Technology Conf., Ottawa, Canada, May 
1998, pp. 2232-2236, space-time coding with OFDM has been 
studied. However, decoding space time codes requires channel 
state information, which is usually difficult to obtain, 
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especially for time-variant channels having dispersive fading. 
Said paper assumes ideal channel state information. 

Another article, also studying space time coding in OFDM having 
assumed known channel state is A. Van Zelst et al. "Space 
5 division multiplexing <SDM) for OFDM systems" presented at VTC 
2000. 

In the article "On channel estimation in OFDMA systems" by Jan- 
Jaap Van de Beek et al. in proceedings IEEE Vehicular Technology 
Conf. Vol. 2, July 1995, page 815-819, methods are disclosed for 

10 channel estimation in OFDM with one transmitter antenna. Minimum 
Mean Square Error (MMSE) and Least Square (LS) estimators with 
reduced complexity are developed that exploit the assumption of 
finite length impulse response. They also address the issue of 
spectral leakage for a discrete time channel impulse response in 

15 a channel having continuous path delays. Considering the aspect 
of single transmit antenna channel estimation, the drawback of 
the proposed methods are that the algorithmic complexity could 
be reduced even further with preserved performance. 

For OFDM systems using space-time coding, two or more different 
20 signals are transmitted from at least two different antennas 
simultaneously. The received signal is the superposition of 
these signals, usually having the same average power. If the 
channel parameters corresponding to each transmitter and 
receiver antenna pair, in a system comprising at least two 
25 transmitting antenna and at least one receiving antennas are 
estimated by the approach developed previously in the article by 
Y. Li et al. "Robust channel estimation for OFDM systems with 
rapid dispersive fading channels IEEE Trans. Commun, Vol. 46, 
pp 902-915, July 1998, the signals from other transmitter 
30 antenna (s) will create interference. The signal-to-interference 
ratio will always be very poor, and the MSE (Minimum Square 
Error) of the estimation will therefore be very large. Hence, 
novel-parameter estimation approaches are desired for 
transmitter diversity using space- time coding. 
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A simple channel estimation for multiple transmit antennas is 
presented in the article by N. Seshadri et al. "Two signaling 
schemes for improving the error performance of frequency- 
selective division duplex transmission systems using transmitter 
5 antenna diversity" in Int. Journal of wireless information 
networks, Vol. 1, No. 1, 1994. The authors propose a method for 
channel learning where one pilot sequence is multiplexed in time 
and over at least two antennas. Hence the pilot is sent on 
antenna 1 first, antenna 2 second and so on. The shortage of 
10 this method is that resources are occupied of sending training 
sequences instead of data. 

In the article by Ye (Geoffrey) Li et al., "Channel estimation 
for OFDM systems with transmitter diversity in mobile wireless 
channels", IEEE Journal on selected areas in communications, 

15 Vol. 17, No. 3, March 1999, the parameter estimation for OFDM 
systems having transmitter diversity is disclosed in detail. In 
said paper transmitter diversity using space-time coding for 
OFDM systems is studied. The channel parameter estimation 
approach is developed. The channel parameter estimation is 

20 crucial for the -decoding of space- time codes, and the MSE bound 
for these estimation approaches is derived. Therefore, OFDM 
systems having transmitter diversity using space-time coding can 
be used for highly efficient data transmission over mobile 
wireless channels. However, the deficiency of the method 

25 disclosed in the above article is that the proposed estimation 
technique is very complex. The structure of the transmitter and 
associated signals is not discussed at all . 

SUMMARY OF THE INVENTION 

A problem with some existing methods for channel estimation in 
30 space time coded systems is that the existing methods consume a 
lot of bandwith since a training symbol is sent for each 
transmit antenna. 
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A further problem with at least one state of the art method 
using concurrent channel training sequences transmissions is 
that a preferred estimation technique is very complex both from 
algorithmic and hardware implementation point of view. 

5 The present invention relates to a method for providing 
multiple channel estimation in the field of space time coding 
in a radio communication system. Space time coding includes 
transmitter diversity, space (-time) multiplexing and other 
complex use of signal coding in time and space. 

10 The inventive solution of the problems is to send pilot signals 
concurrently instead of consecutively in order to increase 
bandwith efficiency during a multitude of channel transfer 
functions in OFDM and doing so, guaranteeing non- interfering 
channel., estimates at low processing cost for both transmitter 

15 and receiver side, in relation to prior art. 

A purpose of the invention is to provide a channel estimation 
technique, which is very bandwith efficient. 

A further purpose of the invention is to provide a channel 
estimation method having a low delay and being very processing 
20 efficient. 

A yet further purpose of the invention is to ensure orthogonal 
and independent channel estimate for up to N transmit antennas, 
the number N being the number of sample values used in the 
receiver and the N channels being frequency flat. 

25 A further purpose of the invention is to reduce hardware 
complexity when providing channel estimation in OFDM- systems. 

A yet further purpose of the invention is to attenuate unwanted 
noise. 

Another purpose of the invention is to be able to provide a low 
30 complexity and robust channel estimation method for a single 
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transmitter antenna case in addition to the more general 
multiple antenna scenario. 

An advantage of the invention is that it is bandwith efficient, 
since merely one synchronization symbol time slot is required 
5 for multiple channel transfer estimations. 

An advantage of the invention is that it is latency and 
processing efficient, since the major parts of the processing 
of the independent channel estimates, per antenna, are 
performed concurrently. 

10 A yet further advantage of the invention is that orthogonality 
is implicitly ensured for up to N transmit antennas by using 
the coding by N times (Inverse) Fast Fourier Transforms 
provided that no channel delay spread exists. This is always 
true disregarding actual content in the training sequence P(K) , 

15 which therefore may be optimized for any other purpose, e.g. 
low peak to average ratio. 

A further advantage of the invention is that hardware 
complexity is very low, since the major part of the processing 
relies on hardware functions inherent for a modem used in a 
20 OFDM system. 

A yet further advantage of the invention is that the channel 
estimates are very exact due to the stuffing with zeros where 
noise is prevalent that acts as a noise attenuation mechanism. 

The term * comprises /comprising" when used in this specification 
25 is taken to specify the presence of stated features, integers, 
steps or components but does not preclude the presence or 
addition of one or more other features, integers, steps, 
components or groups thereof. 

Further scope of applicability of the present invention will 
30 become apparent from the detailed description given 
hereinafter. However, it should be understood that the detailed 
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description and specific examples, while indicating preferred 
embodiments of the invention, are given by way of illustration 
only, since various changes and modifications within the spirit 
and scope of the invention will become apparent to those 
5 skilled in the art from this detailed description. 



BRIEF DESCRIPTION OF THE DRAWINGS 



The invention will now be described in more detail with 
reference to preferred exemplifying embodiments thereof and 
also with reference to the accompanying drawings, in which: 



10 Figure 1 is a schematic view illustrating pilot signals 
being sent separately typical for a state-of-the- 
art method. 

Figure 2a is a schematic view illustrating pilot signals being 
sent concurrently. 
15 Figure 2b is a schematic view illustrating pilot signals being 

sent concurrently from different nodes. 
Figure 3 is a schematic view illustrating the usage of cyclic 

rotation according to the invention. 
Figure 4 is a diagram illustrating a composite channel 
20 impulse response estimate provided with some added 

Gaussian noise. 

Figure 5 is a diagram illustrating further method according 

to the invention. 
Figure 6 is a schematic view illustrating the process at the 
25 receiving antenna. 



DETAILED DESCRIPTION 

In Fig. 1 is depicted a communication system, wherein training 
sequences are transmitted sequentially in an OFDM system 
comprising several transmitting antennas. The training 
3 0 sequences are transmitted in such a way that they do not 
overlap each other over time. The system comprises a 
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transmitter 110 provided with antennas xl, x2, x3 . In figure 1, 
though only three antennas xl, x2, x3 are shown, a man skilled 
in the art understands that more than three antennas may be 
used. The system also comprises at least one receiver 130 
5 provided with at least one receiving antenna yl, y2 . Fig. 1 
illustrates a state-of-the-art approach of performing multiple 
channel estimates in OFDM. Several training symbols pi, p2, p3 
are sent one by one after each other from the different 
antennas xl, x2, x3 to a receiving antenna yl on channels 150, 
10 151, 152. 

In Fig. 2a is illustrated a system according to the invention 
as well as one state of the art method. The system provides a 
multitude of multiple training sequences that are transmitted 
at the same time. The OFDM system in Fig. 2a comprises a 

15 transmitter 210 provided with antennas xl, x2, x3 and at least 
one receiver 230, e.g. a wireless terminal, provided with at 
least one antenna yl . Training sequences pi, p2, p3 are sent 
concurrently from the antennas xl, x2, x3 to the receiving 
antenna yl on channels 250, 251, 252. In contrast to the state 

20 of the art method, the invention uses training sequences pi, 
p2, p3 having a mathematical relationship with each other 
enabling distinguishing estimates in an optimal and simple 
manner. 

In Fig. 2b is illustrated a system according to the invention. 
25 The system provides a multitude of multiple training sequences 

that are transmitted at the same time. The OFDM system in Fig. 

2b comprises a transmitter 610 provided with antennas xl, x2, 

and a transmitter 611 provided with an antenna x3, and at least 

one receiver 230, e.g. a wireless terminal, provided with at 
30 least one antenna yl, y2 . Training sequences pi, p2, p3 are 

sent concurrently from the antennas xl, x2, x3 to the receiving 

antenna yl on channels 250, 251, 252. 

Fig. 3 illustrates transmission of pilot signals according to 
the invention. A known channel estimate sequence P(k) in the 
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frequency domain is provided to a block 308. In the block 308, 
the sequence P(k) is Inverse Fast Fourier Transformed to a 
sequence p(n). The sequence p(n) is fed to a block 360, to a 
block 381 and to a block 382. In the block 360, a cyclic prefix 
5 CP is inserted to precede the sequence p (n) . In a f urther 
embodiment, a cyclic suffix could be used. The cyclic prefix CP 
mitigates intersymbol interference (ISI) effects. Then the 
sequence p(n) is provided to a Digital-to-Analog (D/A) 
converter 370 where it is converted to an analog signal. Then 

10 the D/A converted sequence p(n) is provided to a first antenna 
xl. In the block 381, the sequence p(n) is cyclically rotated 
by a predetermined step comprising a predetermined number n' of 
positions in the sequence and thereby transformed to a sequence 
p(n-n'). Thereby, the sequence p(n-n') is provided to a block 

15 361. In block 361, a cyclic prefix CP is inserted to precede 
the sequence p(n-n'). Thereafter, the sequence p(n-n') is 
provided to a Digital-to-Analog converter 371 where it is 
converted from a digital to an analog shape. Then the D/A 
converted sequence p(n-n') is provided to a second antenna x2 . 

20 The two sequences are then transmitted concurrently from the 
antennas xl and x2 to a wireless terminal 330 provided with at 
least one . antenna yl, y2. In an embodiment comprising three 
antennas, the training sequence p(n) is in the block 382 
cyclically rotated by a second predetermined step comprising 

25 (n"-n') positions in the sequence and thus transformed to a 
sequence p(n''). The sequence p(n-n'') is thus provided to a 
block 362. In block 362, a cyclic prefix CP is inserted to 
precede the sequence p(n-n''). Thereafter, the cyclically 
extended sequence p(n-n'') is provided to a Digital-to-Analog 

30 converter (D/A) 372 where it is converted from a digital to an 
analogue shape. Then the D/A converted sequence p(n-n'') is 
provided to a third antenna x3. The three sequences are then 
transmitted concurrently from the at least one antenna xl, x2 
and x3 to the wireless terminal 330 provided with the antenna 

35 yl. A man skilled in the art understands that the inventive 
system may comprise more than three antennas. 
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The received discrete time signal for two transmit antennas and 
one receive antenna is: 

ri(n)=p(n) *h n (n)+p(n-n' ) *h 2 i (n) +vi (n) , where 

hn and h 2 i are channel impulse responses, and 

5 vi is a noise source associated with the receive antenna. 

A result similar to a cyclic convolution is obtained when an N 
point portion is cut out from an OFDM symbol/ i.e. the removal 
of the cyclic prefix, after the last signal representing the 
previous symbol has arrived. The integer N defines the number 
10 of positions in a Fourier Transformation. 

Then, in the discrete frequency domain, the sequence r x (n) is 
by means of a Fast Fourier Transform transformed to 

R 1 (k)=P(k). Hu(k)+P(k) • e (i " a * n " k}/ " - H 2 i(k)+ Vi(k) 

The sequence Ri(k) can now be divided by the known training 
15 sequence P(k) . The cyclic delay n' dependent phase term remains 
and can mathematically be associated with H 2 i(k) . A subsequent 
Inverse Fast Fourier transform returns a composite channel 
impulse response estimate, 

hconp (n) =hn (n) +h 2 i (n-n' ) +Vi (n) 

20 Now, a commonly used constraint in OFDM is that the cyclic 
prefix length should be slightly longer than the longest 
duration of a channel impulse, i.e. delay spread. In addition, 
the OFDM symbol duration is often chosen in order to enable the 
cyclic prefix to add maximally 20-25% overhead, but preferably 

25 less than that as both energy and efficiency is consumed. 
Provided that the number of cyclic positions shift n' between 
the two antenna paths is greater than the length of the cyclic 
prefix CP, the channel impulse responses hn and h 2 i are 
guaranteed to . be extracted individually from a composite 

30 channel impulse response estimate hc^p.. An example of the 
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magnitude (or real or imaginary part) in dB of a composite 
channel impulse response estimate Ihcompl having some added 
Gaussian noise is indicated in Fig. 4i In the example 
illustrated in Fig. 4 h u and h 21 are 4 and 3 point long 
5 respectively, and N=32 . 

The individual channel impulse responses are easily extracted 
from the received signal since the cyclic shift of n' positions 
is known, and the extraction can be made as illustrated by 
examples shown in Figs. 4 and 5. The cyclic shift of n' 

10 positions for the individual channel impulse responses must be 
removed prior they are used for equalization of the channels. 
Then, the channel impulse response corresponding to h2i is 
cyclically shifted n' positions in the opposite direction 
compared to the direction in which the associated pilot 

15 sequence originally was shifted, as this unfolds the phase 
wrapping on the corresponding channel transfer function, else 
being incurred by the cyclic shift. A cyclic shift is performed 
in a vector having N positions, the last position, or the two 
last positions being put at the beginning of the vector, and 

20 the other positions being shifted to the right, clockwise. In a 
further embodiment the positions are shifted counterclockwise. 
For each individual channel impulse response, positions being 
considered irrelevant or having non- significant energy are 
replaced with zero values, alternatively, each individual 

25 channel impulse response being more smoothly attenuated by a 
smother windowing function/ see Figs. 4 and 5. A man skilled in 
the art understands that the replacement of only noise 
dependent parts of the channel impulse responses with zero 
results in a significant reduction of noise and therefore 

30 provides for excellent estimation. For each channel estimate, 
noise may be reduced with roughly 10-logio (N/delay spread) dB. 
The delay spread indicates the extent of the energy dispersion 
from different multipath components, i.e. the greater the delay 
spread, the later arrives the last multipath component . 
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Two methods for separating two channel impulses responses will 
now be described with reference to Figs. 4 and 5. The channel 
impulse response has about the same length as the cyclic 
prefix. 

The first method is illustrated by the diagram in fig. 4. 
According to the first method, a plurality of fixed ranges is 
defined in the discrete time domain. The number of fixed ranges 
is equal to the number of transmit antennas and are upper 
limited to floor (N/CP duration) . Floor is a mathematical 
function choosing the closest lower integer value of the 
argument. From each of said ranges, one channel impulse 
response is extracted. For example, in range 1, a number of 
positions have amplitudes in the composite channel impulse 
sequence significantly stronger than other positions. The 
strong channel impulse responses are of interest for the 
inventive method. 

As the cyclic shifts are well known to both the transmitter and 
the receiver, the receiver knows when to expect the start of 
the ranges . As the receiver also knows the length of the cyclic 
20 prefix, it also knows the end point of each range. The span of 
the range is at least the length of the cyclic prefix having 
additional extension to protect from leakage. It is assumed 
that the timing synchronization for . the Fast Fourier 
Transformation window start has been performed at the preceding 
25 timing synchronization step. However, small errors in timing 
are acceptable, due to the fact that said errors are transposed 
to a small time shifts in the composite channel impulse 
response and thus each channel impulse response is kept within 
each individual range. Very large timing synchronization errors 
30 may require larger ranges in order to accommodate protection 
for such imperfections in preceding steps. 

It should be noted that the leakage and significant timing 
errors may rotate into the end of the composite channel impulse 
response due to cyclic properties of the Fast Fourier 
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Transformation. The leakage explains why range 1 extends 
cyclically into the end of the composite channel impulse 
response. When creating the channel impulse response 
corresponding to range 1 in the composite channel impulse 
5 response, positions of range 1 are copied, whereas the 
positions outside the range 1 are set to 0. The positions 
outside range 1 have insignificant contributions to the channel 
impulse response hn and may therefore be set to zero. 

For channel impulse response within range 2, which are the 
10 results of a rotation of the training sequence p by n' 
positions, the extraction follows the same procedure as for the 
range 1. However, the resulting channel impulse response 
sequence for range 2 is cyclically shifted back n' positions in 
such a way that most of the energy is gathered at the beginning 
15 of the channel impulse response. 

For an embodiment using more than two antennas at the 
transmitter end, the same procedure is performed, but having 
different range values and subsequent cyclic rotation. Those 
values depend on how the cyclic shift is performed for the 
20 transmission of respective cyclic rotated training sequence. 

A further method is illustrated by the diagram of Fig 5. This 
method comprises an additional step. This method employs a more 
intelligent and adaptive mechanism. For each range, a first 
position having the maximum magnitude within the fixed range is 

25 determined. Thereafter, second and third positions are 
determined, being the furthermost positions within the range on 
the respective side and away ' from the first position also 
having a predetermined lower threshold with respect to the 
first position. The adaptively chosen subset from the fixed 

30 range is copied, whereas the remaining positions are set to 
zero. This step is performed for each fixed range. The figure 
shows the case for the fixed range 2. A cyclic rotation is 
performed in the same way as the embodiment disclosed in 
connection with Fig. 4. It should be noted that if all 
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magnitudes are within the predefined magnitude difference, the 
entire range will stay unaltered. 

A man skilled in the art understands that the above disclosed 
scheme may be used for a greater number of transmit antennas. 
5 However,, at most the cyclic prefix CP, floor (N/CP duration), 
undisturbed channel transfer function, i.e. transmission over 
channels having a minimum of interference, can be guaranteed 
per receiving antenna. In addition, some guard positions must 
be left between the individual channel impulse responses due to 
10 leakage from non-whole point delays. 

In a further embodiment of the invention, additionally 
windowing is used in order to minimize the effect of Fast 
Fourier Transform leakage when path delays are intermediately 
placed between sampling points. Leakage is an inherent 

15 characteristic of the (I)DFT process when the real, i.e. 
continuous, world is sampled and results in some interference 
leakage to neighboring sample points. One such real world 
effect is e.g. that path delays may be somewhere intermediately 
placed between the sampling instances. In one embodiment of the 

20 windowing, a window function is applied in the frequency domain 
after the division by the sequence P(k), but prior to the IFFT 
to the composite channel impulse response and subsequent 
channel impulse response partitioning. After the partitioning 
each channel estimate should be compensated for the effects of 

25 the performed windowing by performing inverse windowing. 

A well known window function is the Hanning window also called 
Raised cosine window. In a preferred embodiment, the Hanning 
windowing in the frequency domain is replaced by a 
corresponding time domain filter operation on the composite 
30 channel impulse response according to hconp_ W indowed (n) =h CO inp(n) /2- 
hconp(n-l)/4-hc OI np(n+l)/4, where the delay corresponds to cyclic 
shifts of hcompfn). This signal processing trick allows low 
complexity implementation of the frequency domain Hanning 
window operation. 
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In an embodiment of the invention comprising four transmitting 
antennas the composite channel impulse response estimate 
becomes 

Hcoinp (n) =hu (n-ni ' ) +h 2X (n-n 2 ' ) +h 3 i (n-n 3 ' ) +h 4 i (n-n* ' ) +wi (n) 

5 n defines the cyclic rotation, wherein ni'=£n 2 '*n 3 '?tn4' , and 
preferably ni '= 0, n 2 '> ni'+ CP_duration, n 3 '> n 2 ' + CP_duration, 
n 4 '> n 3 ' + CP — duration and n 4 '< N-CP_duration, 

It should be noted that only two Fast Fourier Transformations 
are required in block 666, fig. 6, i.e. one Inverse and one 
10 ordinary Fast Fourier Transf ormation. In a further embodiment, 
the individual extracted channel impulse responses are Fast 
Fourier Transformed to their respective frequency domain 
representation Hxy(k) under the assumption that block 634 in 
fig. 6 operates with channel transfer functions. 

15 Practically a full Fast Fourier Transf ormation is used. 
However, in a further embodiment, for each individual channel 
impulse response, the last Fast Fourier Transformation is of a 
reduced type, since the major content is zeros. A traditional 
(Inverse) Fast Fourier transformation can remove some 

20 operations, when part of the input signals are known to 
comprise zeros, as is well known for a man skilled in the art. 
Alternatively the frequency response may be generated directly 
from the DFT definition when the number on non-zero elements in 
hcomp (n) are few. 

25 Moreover, if the receiver exploit multiple receiving antennas, 
as is indicated in Fig. 3 the inventive process can be 
performed independently for each antenna. 

Now, a receiver structure for OFDM will be disclosed with 
reference to Fig. 6. In Fig. 6, only one antenna is 
30 illustrated. However, a man skilled in the art understands that 
multiple antennas may be used. Signals are received at the 
antennas yl, y2... yJ, where J is an integer ranging from 1 and 
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upwards limited only by practical concerns such as space, power 
consumption, cost etc. In Fig. 6 only the first antenna yl is 
shown. The signal received at the antenna yl is provided to a 
block 612, where it is amplified to an appropriate level for 
5 reception by means of an amplifying stage, often with Automatic 
Gain Control (AGC) . Then, the amplified signal is provided to a 
block 614, being an A/D converter, where the signal is 
converted from an analog shape to a digital shape and samples 
are taken of the received amplified signal at successive times. 

10 Then, the digital signal is provided to a block 616, where 
frequency and OFDM timing synchronization is performed. The 
timing synchronization is often based on particular training 
symbols for this specific purpose, but other methods known for 
the man skilled in the art may be used. When the timing 

15 synchronization has been performed, the signal is provided to a 
block 618, where the cyclic prefix is removed. A discrete and 
complex valued sequence r P ,i(n) comprising training information 
is provided from the block 618 to a channel estimation block 
650, which will be disclosed in more detail below. A signal 

20 r s ,i(n) comprising modulated data from multiple transmit 
antennas and affected by the channel is provided from the block 
618 to a block 620, where a Fast Fourier transformation is 
performed. The block 620 provides a signal R s ,i(k) for each 
receiving antenna yl, y2 ... yj per OFDM symbol interval. The 

25 signal R Bl i(k) is provided to a demodulation/ decoding block 634, 
where channel equalization, de-interleaving, FEC (forward error 
correction) decoding and de-scrambling is provided. 

In order to decode the data that will be sent, the channel need 
to be estimated. The channel estimation is made in the block 

30 650. The signal r P/X (n) from block 618 is provided to a block 
660, where a Fast Fourier Transformation is performed resulting 
in a signal R P ,i(k). The signal R P/ i(k) is provided to a block 
662, where the signal R P/i (k) is divided by a training symbol 
P(k). The result from block 662 is provided to a block 664, 

35 where an Inverse Fast Fourier Transformation is performed. When 



WO 02/45329 



PCT/SE01/02292 



22 



the training symbol sequence arrives, the time domain signal is 
switched into the channel estimation stage, where a composite 
channel impulse response h'(n) is extracted. The composite 
channel impulse response h'(n) is then fed from the block 664 
5 to a block 666. In block 666, individual channel impulse 
responses for this particular reception antenna yl and all 
transmit antennas are then extracted and fed to the 
decoding/demodulation block 634. The decoding/demodulation 
block 634 uses the channel estimates when the signals Rs,i(k), 
10 R s ,2(k),... / R s ,j(k) are processed. The signal processed in block 
634 is then output 630 to higher layers. 

Depending of the operation with respect to time or frequency in 
the decoding/demodulation block 634, it should be noted that 
often a frequency domain representation of the channel may be 

.15 more suitable than the time domain oriented channel impulse 
response. Block 634 may for example use a frequency domain 
equalizer prior Forward Error Correction decoding. 
Alternatively, block 634 may incorporate frequency domain . 
channel state information in a Viterbi decoder's path metric 

20 when convolutional coding is employed. Hence, additional FFTs 
are required for this purpose. Furthermore, said operation is 
performed in parallel for potential multiple receiver antennas. 
It should also be noted that the Fast Fourier Transform 
function in block 620, 660 and additional blocks having an FFT 

25 mechanism of the same size may use the same hardware. It should 
be noted that the disclosed invention in addition to the 
channel estimate capability also enables single channel 
estimates from a single transmitter antenna. The Inverse Fast 
Fourier function may also use the same hardware with small 

30 alterations of involved multiplication factors. 
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CLAIMS 

1. A method in a wireless communication system comprising at 
least one transmitter (210) provided with at least one 
antenna (xl, x2, x3) and at least one receiving unit (330) 

5 provided with at least one antenna (yl, y2) and wherein 

training sequences (pi, p2, p3) are transmitted from the at 
least one antenna (xl, x2, x3) of the at least one 
transmitter (210) to the at least one antenna (yl, y2) of 
the at least one mobile unit, characterized in that 

10 - first, prior to the transmission, a training sequence P(k) is 
Inverse Discrete Fourier Transformed to a sequence p(n); 

- second, for each antenna branch the Inverse Discrete Fourier 
Transformed sequence p(n) is cyclically rotated by a number 
of predetermined steps (m, n 2 ) being different for each 

15 antenna branch; 

- third, the cyclically rotated training sequences p(n-ni), 
p(n-n 2 ) are transmitted concurrently from different antennas 
to the receiving unit; and 

- fourth, at the receiving unit, the received sequences (s(n- 
20 ni) , s(n-n 2 )) being a superposition of transmitted training 

sequences, each individually affected by the propagation 
medium, are used to provide channel impulse response 
estimates for the transmission from respective antenna. 

2. A method in a wireless communication system comprising at 
25 least one transmitter (210) provided with at least two 

antennas (xl, x2, x3) and at least one receiving unit (330) 
provided with at least one antenna (yl, y2) and wherein 
training sequences (pi, p2 , p3) are transmitted from the at 
least two antennas (xl, x2, x3) of the at least one 
30 transmitter (210) to the at least one antenna (yl, y2) of 

the at least one mobile unit, characterized in that 
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- first, prior to the transmission, a training sequence P(k) is 
Inverse Discrete Fourier Transformed to a sequence p(n) ; 

- second, for each antenna branch the Inverse Discrete Fourier 
Transformed sequence p(n) is cyclically rotated by a number 
of predetermined steps (ni, n 2 ) being different for each 
antenna branch; 

- third, the cyclically rotated training sequences p(n-ni), 
p(n-n2> are transmitted concurrently from different antennas 
to the receiving unit; and 

- fourth, at the receiving unit, the received sequences (s(n- 
ni) , s(n-n 2 )) being a superposition of transmitted training 
sequences, each individually affected by the propagation 
medium, are used to provide channel impulse response 
estimates for the transmission from respective antenna. 

3. A method according to claim 2, characterized in that instead 
of a Inverse Discrete Fourier Transform, an Inverse Fast 
Fourier Transform is performed. 

4. A method according to claim 2 or 3, characterized in that 

- a cyclic extension having a predetermined length is added to 
each sequence prior the transmission, the cyclic extension 
being greater than the delay spread. 

5. A method according to claim 2 or 3, characterized in that 
the distance between each step is greater than the delay 
spread . 

6. A method according to claim 2, characterized in that at the 
receiving unit the received sequence 

- in a first step, is Discrete Fourier Transformed and divided 
by the training sequence P(k), 
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- in a second step, the result from the first step is Inverse 
Discrete Fourier Transformed resulting in a sequence having 
distinctly separated regions in the time domain, the 
separated regions containing the respective channel impulse 

5 response estimates. 

7. A method according to claim 3, characterized in that at the 
receiving unit the received sequence 

- in a first step, is Fast Fourier Transformed and divided by 
the training sequence P(k), 

10 - in a second step, the result from, the first step is Inverse 
. Fast Fourier Transformed resulting in a sequence having 
distinctly separated regions in the time domain* the 
separated regions containing the respective channel impulse 
response estimates. 

15 8. A method according to claim 6 or 7, characterized in that 

- fixed predetermined ranges are selected in the discrete time 
domain, 

- each range comprising one and only one of the above defined 
regions, 

20 - one channel impulse response is selected from each of said 
ranges , 

9. A method according to claim 8, characterized in that 

- in each range, only the strongest positions are used, and the 
others are set to zero, and, for each selected channel 

25 impulse response, positions outside the range are replaced by 

zeros , 

10. A method according to claim 8 or 9, characterized in that 
each resulting channel impulse response is converted to the 
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frequency domain by a DFT or FFT operation depending on 
operating methods of primarily equalizer and FEC decoder. 

11. A method according to claim 6 or 7, characterized in that a 
window function is applied prior to the second step, wherein 

5 leakage inherent in the preceding transformation in the 

first step is reduced. 

12. A method according to claim 6 or 7, characterized in that a 
filter function is applied after the second step, wherein 
leakage inherent in the preceding transformation in the 

10 first step is reduced. 

13. A method according to claim 11 or 12, characterized in that 
said window function is a Harming window or said filter 

< Inversion is an IDFT transformed Hanning window. 

14. A method according to any of claims 8, 9, 11 or 12, 
15 characterized in that an inverse impulse response 

corresponding to the window function is applied after the 
selection of said channel impulse response, wherein the 
phase and amplitude values are compensated due to the result 
from the preceding window. 

20 15. A method according to any of claims 8-12, characterized in 
that an inverse function is applied after the conversion to 
the frequency domain by a DFT or FFT operation, wherein the 
phase and amplitude values are compensated due to the result 
from the preceding window. 

25 16. A wireless communication system comprising at least one 
transmitter (210) provided with at least two antennas (xl, 
x2, x3) and at least one receiving unit (330) provided with 
at least one antenna (yl, y2) and wherein training sequences 
(pi, p2, p3) are transmitted from the at least two 

30 antennas (xl, x2 , x3) of the at least one transmitter (210) 

to the at least one antenna (yl, y2) of the at least one 
receiver unit, characterized in that 
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- first, prior to the transmission, a training sequence P(k) is 
Inverse Discrete Fourier Transformed to a sequence p(n); 

- second, for each antenna branch the Inverse Discrete Fourier 
Transformed sequence p(n) is cyclically rotated by a number 

5 of predetermined steps {n lt n 2 ) / being different for each 

antenna branch; 

- third, the cyclically rotated training sequences p(n-ni), 
p(n-n 2 ) are transmitted concurrently from different antennas 
to the receiving unit; and 

10 - fourth, at the receiving unit, the received sequences (s(n- 
ni), s(n-n 2 )) being a superposition of transmitter training 
sequences each individually affected by the propagation 
medium, are used to provide channel impulse response 
estimates for the transmission from respective antenna. 

15 17. A system according to claim 16, characterized in that 
instead of a Inverse Discrete Fourier Transform, an Inverse 
Fast Fourier Transform is performed. 

18. A system according to claim 16 or 17, characterized in that 

- a cyclic extension having a predetermined length is added to 
20 each sequence prior the transmission, the cyclic extension being 
greater than the delay spread. 

19. A system according to claim 16 or 17, characterized in 
that the distance between each step is greater than the delay 
spread. 

25 20. A system according to claim 16, characterized in that at the 
receiving unit the received sequence 

- in a first step, is Discrete Fourier Transformed and divided 
by the training sequence P(k), 
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- in a second step, the result from the first step is Inverse 
Discrete Fourier Transformed resulting in a sequence having 
distinctly separated regions in the time domain, the 
separated regions containing the respective channel impulse 

5 response estimates. 

21. A system according to claim 17, characterized in that at the 
receiving unit the received sequence 

- in a first step, is Fast Fourier Transformed and divided by 
the training sequence P(k) , 

10 - in a second step, the result from the first step is Inverse 
Fast Fourier Transformed resulting in a sequence having 
distinctly separated regions in the time domain, the 
separated regions containing the respective channel impulse 
respons.e estimates. 

15 22. A system according to claim 20 or 21, characterized in that 

- fixed predetermined ranges are selected in the discrete time 
domain, 

- each range comprising one and only one of the above defined 
regions, 

20 - one channel impulse response is selected from each of said 
ranges , 

23. A system according to claim 22, characterized in that 

- in each range, only the strongest positions are used, and the 
others are set to zero for each selected channel impulse 

25 response positions outside the range are replaced by zeros. 

24. An arrangement in a wireless communication system comprising 
at least one transmitter (210) provided with at least two 
antennas (xl, x2, x3) and at least one receiving unit (330) 
provided with at least one antenna (yl, y2) and wherein 
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training sequences (pi, p2, p3) are transmitted from the at 
least two antennas (xl, x2, x3) of the at least one 
transmitter (210) to the at least one antenna (yl, y2) of the 
at least one receiver unit, characterized by 

5 - means for performing, prior to the transmission, an Inverse 
Discrete Fourier Transformation on a training sequence P(k) 
to produce a sequence p (n) ; 

- means for performing, for each antenna branch, a cyclic 
rotation by a predetermined number of steps (n x , n 2 ) of the 

10 Inverse Discrete Fourier Transformed sequence, said steps 

being different for each antenna branch; 

- means for transmitting concurrently the cyclically rotated 
training sequences p(n-ni) , p(n-n 2 ) from different antennas to 
the receiving unit; and 

15 - means for using, at the receiving unit, the received 
sequences (s(n-ni), s(n-n 2 )), being a superposition of 
transmitted training sequences each individually affected by 
the propagation medium, to provide channel impulse response 
estimates for the transmission from respective antenna. 

20 2 5. An arrangement according to claim 23, characterized in that 
instead of a Inverse Discrete Fourier Transform, an Inverse 
Fast Fourier Transform is performed. 

26. An arrangement according to claim 24 or 25, characterized 
by means for adding a cyclic extension having a 

25 predetermined length to each sequence prior the 

transmission, the cyclic extension being greater than the 
delay spread. 

27. An arrangement according to claim 24 or 25, characterized 
in that the distance between each step is greater than the 

30 delay spread. 



WO 02/45329 PCT/SEO 1/02292 

30 



28. An arrangement according to claim 24, characterized in that 
the receiving unit comprises means for 

- perfo:nrning, in a first step, a Discrete Fourier Transform and 
a division of the received sequence by the training sequence 

5 P(k), 

- performing, in a second step, an Inverse Discrete Fourier 
Transform of the result from the first step resulting in a 
sequence having distinctly separated regions in the time 
domain, the separated regions containing the respective 

10 channel impulse response estimates. 

29. An arrangement according to claim 25, characterized in that 
the receiving unit comprises means for 

- performing, in a first step, a Fast Fourier Transform and a 
division of the received sequence by the training sequence 

15 P(k), 

- performing, in a second step, an Inverse Fast Fourier 
Transform of the result of the first step resulting in a 
sequence having distinctly separated regions in the time 
domain, the separated regions containing the respective 

20 channel impulse response estimates. 

30. An arrangement according to claim 28 or 29, characterized in 
that the receiving unit comprises 

- means for selecting fixed predetermined ranges in the 
discrete time domain, each range comprising one and only one 

25 of the above defined regions, and 

- means for selecting one channel impulse response from each of 
said ranges, 

31. An arrangement according to claim 30, characterized in that 



WO 02/45329 PCT/SEO 1/02292 

31 

- in each range, only the strongest positions are used, and the 
others are set to zero for each selected channel impulse 
response position outside the range are replaced by zero. 
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